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Effect of Crystal Orientation on Fatigue Failure of Single Crystal Nickel Base Turbine 

Blade Superalloys 


1.0 Introduction 

High Cycle Fatigue (HCF) induced failures in aircraft gas-turbine engines is a pervasive problem 
affecting a wide range of components and materials. HCF is currently the primary cause of 
component failures in gas turbine aircraft engines \ Turbine blades in high performance aircraft 
and rocket engines are increasingly being made of single crystal nickel superalloys. Single- 
crystal Nickel-base superalloys were developed to provide superior creep, stress rupture, melt 
resistance and thermomechanical fatigue capabilities over polycrystalline alloys previously used 
in the production of turbine blades and vanes. Currently the most widely used single crystal 
turbine blade superalloys are PWA 1480/1493 and PWA 1484. These alloys play an important 
role in commercial, military and space propulsion systems. PWA1493, identical to PWA 1480, 
but with tighter chemical constituent control, is used in the NASA SSME alternate turbopump, a 
liquid hydrogen fueled rocket engine. 

Single-crystal materials differ significantly from polycrystalline alloys in that they have highly 
orthotropic properties making the position of the crystal lattice relative to the part geometry a 
significant factor in the overall analysis 2 . The modified Goodman approach currently used for 
component design does not address important factors that affect HCF such as fretting and/or 
galling surface damage, and interaction with LCF \ Rocket engine service presents another set of 
requirements that shifts emphasis to low temperature fatigue and fracture capability with 
particular attention given to thermal, cryogenic and high pressure hydrogen gas exposure \ To 
address HCF induced component failures, the gas turbine industry, NASA, the US Air Force, and 
the US Navy have made significant efforts in understanding fatigue in single-crystal turbine 
blade superalloys. Understanding fatigue initiation, threshold, and Region II fatigue crack growth 
are of primary importance and there is great need for improvements in fracture mechanics 
properties of turbine blade alloys. While a large amount of data has been collected there 
currently exists no simple method for applying this knowledge toward the design of more robust 
single crystal gas turbine engine components. It is therefore essential to develop failure criteria 
for single crystals, based on available fatigue and fracture test data that will permit a designer to 
utilize the lessons learned. 

Objectives for this paper are motivated by the need for developing failure criteria and fatigue life 
evaluation procedures for high temperature single crystal components, using available fatigue 
data and finite element modeling of turbine blades. Fatigue failure criteria are developed for 
single crystal material by suitably modifying failure criteria for polycrystalline material. The 
proposed criteria are applied for uniaxial LCF test data, to determine the most effective failure 
parameter. A fatigue life equation is developed based on the curve fit of the failure parameter 
with LCF test data. Single crystal turbine blades used in the alternate advanced high-pressure 
fuel turbopump (AHPFTP/AT) are modeled using a large-scale 3D finite element (FE) model 
capable of accounting for material orthotrophy and variation in primary and secondary crystal 
orientation. Using the FE stress analysis results and the fatigue life relations developed, the effect 
of variation of primary and secondary crystal orientations on life is determined, at critical blade 
locations. The most advantageous crystal orientation for a given blade design is determined. 
Results presented demonstrates that control of secondary and primary crystallographic 



Nagaraj K. Arakere and Gregory R. Swanson 


Page 3 


orientation has the potential to optimize blade design by increasing its resistance to fatigue crack 
growth without adding additional weight or cost. 

2.0 Crystal Orientation 

Nickel based single-crystal superalloys are precipitation strengthened cast mono grain 
superalloys based on the Ni-Cr-Al system. The microstructure consists of approximately 60% by 
volume of y* precipitates in a y matrix. The y’ precipitate, based on the intermetallic compound 
N^Al, is the strengthening phase in nickel-base superalloys. The / precipitate is a Face Centered 
Cubic (FCC) structure and composed of the intermetallic compound N^Al. The / precipitate is 
suspended within the y matrix, which is also of FCC structure and comprised of nickel with 
cobalt, chromium, tungsten and tantalum in solution. Single crystal superalloys have highly 
orthotropic material properties that vary significantly with direction relative to the crystal lattice. 
Primary crystallographic orientation of a turbine blade, commonly referred to as a, is defined as 
the relative angle between the airfoil stacking line and the <001> direction, as shown in Fig. 1. 
Current manufacturing capability permits control of a to within 5° of the stacking line. 
Secondary orientation (i defines the angle of the <100> orientation relative to the blade 
geometry. In most turbine blade castings the secondary orientation (3 is neither specified nor 
controlled during the manufacturing process. The P orientation for a given blade casting 
therefore becomes a random variable. Usually, however, the p orientation for each blade is 
recorded after the casting process is complete. 

3.0 Fatigue in Single Crystal Nickel Superalloys 

Slip in metal crystals often occurs on planes of high atomic density in closely packed directions. 
The 4 octahedral planes corresponding to the high-density planes in the FCC crystal are shown in 
Fig. 2. Each octahedral plane has 6 slip directions associated with it. Three of these are termed 
easy-slip or primary slip directions and the other 3 are secondary slip directions. Thus there are 
12 primary and 12 secondary slip directions associated with the 4 octahedral planes 4 . In 
addition, there are 6 possible slip directions in the three cube planes, as shown in Fig. 3. 
Deformation mechanisms operative in PWA 1480/1493 are divided into three-temperature 
regions 5 . In the low temperature regime (26C to 427C, 79F to 800F)) the principal deformation 
mechanism is by (111 )/< 1 1 0> slip, and hence fractures produced at these temperatures exhibit 
(111) facets. Above 427C (800F) thermally activated cube cross slip is observed which is 
manifested by an increasing yield strength up to 87 1C (1600F) and a proportionate increase in 
(111) dislocations that have cross slipped to (001) planes. Thus nickel based FCC single crystal 
superalloys slip primarily on the octahedral and cube planes in specific slip directions. At low 
temperature and stress conditions crystallographic initiation appears to be the most prevalent 
mode. This mode warrants special consideration since this mode of cracking has been observed 
in many turbine blade failures \ The operative deformation mechanism has a strong influence on 
the nature of fracture. As a result of the two phase microstructure present in single crystal nickel 
alloys a complex set of fracture mode exists based on the dislocation motion in the matrix (ft and 
precipitate phase (/). Telesman and Ghosn 6 have observed the transition of fracture mode as a 
function of stress intensity (K) in PWA1480 at room temperature. Deluca and Cowles 7 have 
observed the fracture mode transition that is environmentally dependent. 
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Fatigue life estimation of single crystal blade components represents an important aspect of 
durability assessment. Turbine blade material is subjected to large mean stresses from the 
centrifugal stress field. High frequency alternating fatigue stresses are a function of the vibratory 
characteristics of the blade. Any fatigue failure criteria chosen must have the ability to account 
for high mean stress effects. Towards identifying fatigue failure criteria for single crystal 
material we consider four fatigue failure theories used for polycrystalline material subjected to 
multiaxial states of fatigue stress. Kandil, et al 8 presented a shear and normal strain based 
model, shown in Eq. 1, based on the critical plane approach which postulates that cracks initiate 
and grow on certain planes and that the normal strains to those planes assist in the fatigue crack 
growth process. In Eq. (1) y max is the max shear strain on the critical plane, £„ the normal strain 
on the same plane, S is a constant, and N is the cycles to initiation. 

)W+S£„=/W (1) 

Socie, et al 9 presented a modified version of this theory shown in Eq. (2), to include mean stress 
effects. Here the maximum shear strain amplitude (A# is modified by the normal strain 
amplitude (Ac) and the mean stress normal to the maximum shear strain amplitude (cg to ). 


a y | a g„ 
2 2 



= f(N) 


( 2 ) 


Fatemi and Socie 10 have presented an alternate shear based model for multiaxial mean-stress 
loading that exhibits substantial out-of-phase hardening, shown in Eq. (3). This model indicates 
that no shear direction crack growth occurs if there is no shear alternation. 


. max 

+ — -) = /(AO (3) 

2 c 7 y 

Smith, et al 11 proposed a uniaxial parameter to account for mean stress effects which was 
modified for multiaxial loading, shown in Eq. (4), by Bannantine and Socie 12 . Here the 
maximum principal strain amplitude is modified by the maximum stress in the direction of 
maximum principal strain amplitude that occurs over one cycle. 

A p 

m “) = /(W) (4) 

4.0 Application of Failure Criteria to Uniaxial LCF Test Data 

The polycrystalline failure parameters described by Eqs. (1-4) will be applied for single crystal 
uniaxial strain controlled LCF test data. Transformation of the stress and strain tensors between 
the material and specimen coordinate systems (Fig. 4) is necessary for implementing the failure 
theories outlined. The components of stresses and strains in the (x‘, y‘, z‘) system in terms of the 
(x, y, z) system is given by 13 

M=[ fi'JM; fc’}=[£M (5) 


{a} = \QV W) = \Q\ W) ; {<?} = \Q'eV {e'} = fej M 


(6) 
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Where 


M 


ft a 2 ft la, a., 2a, a, 2a, a, 

ft ft ft Iftft 2 ft B, 2 ftP, 

\q}= 7 ' 72 7 ' 2hh l7n 2727 

Pj i P2Y2 PiYs (P2Y2 + PiYi) (An + A/i) (Pji+PiYi) 

m Yi a 2 Y-fih, (Yi a 2 + Ys a i) (Y& + Y&) (YPi+YM) 

cc x P\ a 2 Pi a A (UiPi + ohPi) A Ps + ccsPJ (a l fi 2 +a 2 P) 

and 

a] a 2 a] a 2 a 2 ayx 2 a 2 a, 

Px PI PI PA PA P2P 

\o]= 7; 72 72 7 771 772 7l7{ 

2 P{Y\ 2 PiYi ZPiY) (P2Y2 + P2Y2) (PxYs + PlYx) (PJ2 + P2Y) 

2 Y& 2 / 2 <z 2 2 y,a, (y 2 a 2 + y 2 a 2 ) (y l a i + y i a l ) (ya 2 + y 2 a,) 
la 2 p 2 2a A ( a 2 P\ -t* a^ P 2 ) ( (x, f5 2 ft cc ) (a^ fi 2 ft oc^p ^ ) 

Table 1 shows the direction cosines between the (x, y, z) and (x\ y’, z’) coordinate axes. The 
transformation matrix [ Q ] is orthogonal and hence [Q\ x = [£?] T = [Q ’]. The generalized 
Hooke’s law for a homogeneous anisotropic body in Cartesian coordinates (x, y, z) is given by 
Eq. (10) 13 . 

fc}=MM (io) 

[fly] is the matrix of 36 elastic coefficients, of which only 21 are independent, since [fly] =[0,,]. 
The elastic properties of FCC crystals exhibit cubic symmetry, also described as cubic syngony. 
Materials with cubic symmetry have only three independent elastic constants designated as the 
elastic modulus, shear modulus, and Poisson ratio 4 and hence [fly] has only 3 independent elastic 
constants, as given below. 
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a n a n a 12 0 0 0 

a n a u a n 0 0 0 

a n a n a n 0 0 0 

0 0 0 0 0 

0 0 0 0 a M 0 

0 0 0 0 0 a u 


( 12 ) 


The elastic constants are 


1 _ 1 

— » «44 ~7 r 

G v 


The elastic constants in the generalized Hooke’s law of an anisotropic body, [fly], vary with the 
direction of the coordinate axes. In the case of an isotropic body the constants are invariant in 
any orthogonal coordinate system. The elastic constant matrix [fl’y] in the (x\ y\ z’) coordinate 
system that relates {£’ } and {<f} is given by the following transformation 13 . 


Wy]= M klfe] = J J ^ 4 a mn Q m iQnj 

m — 1 n = 1 

(ij= 1,2 ,6) 


(14) 


Shear stresses in the 30 slip systems shown in figures 1 and 2 are denoted by t 1 , x 2 ,..., x 30 . The 
shear stresses on the 24 octahedral slip systems are 4 



and the shear stresses on the 6 cube slip systems are 
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I &yt 


(16) 


Shear strains (engineering) on the 30 slip systems are calculated using similar kinematic 
relations. 


As an example problem we consider a uniaxial test specimen loaded in the [213] direction 
(chosen as the x’ axis in Fig. 4) under strain control. The applied strain for the specimen is 1.212 
%. We wish to calculate the stresses and strains in the material coordinate system and the shear 
stresses on the 30 slip systems. The x’ axis is aligned along the [213] direction. The required 
direction cosines are shown in Table 2. The stress-strain relationship in the specimen coordinate 
system is given by 


I44;IM 


(17) 


The [a’jj] matrix is calculated using Eq. (14). 



3.537H - 8 

- 2.644E - 9 

-1.986E - 8 

5.209E - 9 

1.405E - 8 

1.878E - 8 

- 2.644E - 9 

3.975E - 8 

-2.423E - 8 

-5.61E-9 

1.297E - 8 

-2.023E - 8 

- 1.986E - 8 

-2.423E - 8 

5.696E - 8 

4.007E - 10 

-2.703E - 8 

1.445E - 9 

5.209E - 9 

-5.61E - 9 

4.007 E - 10 

7.089E - 8 

2.889E - 9 

2.595E - 8 

1.405E - 8 

1.297E - 8 

-2.703E - 8 

2.889E - 9 

8.838E - 8 

1.042E - 8 

1.878E - 8 

-2.023E - 8 

1.445E - 9 

2.595E - 8 

1.042E - 8 

1.572E - 7 


(18) 


Since c x ' is the only nonzero stress in the specimen coordinate system, we have 


a = 


0.01212 

3.537E-8 


= 342,663 psi 


(19) 


Knowing { cf } we can now calculate { e ' } as 


G x 


'342,663' 


0.01212 

< 


0 


-9.059E-4 

e\ 

-w 

0 


-6.805E-3 

< , 


» z= 1 


r n 

0 


1.785E-3 

K 


0 


4.815E-3 

/J 


0 


6.435E-3 

V. 


( 20 ) 
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The stresses and strains in the material coordinate system can be calculated using Eqs. (6) as 




- 1.43E - 5 ' 


a* 


9.789E + 4 

£ y 


-6.693E-3 


°y 


2.447E + 4 

£ z 


0.011 


<y z 


2.203E + 5 


> = < 


> ^ < 


> — < 


Yyi 


4.676E-3 


7z 


7.342E + 4 

r* 


9.353E-3 




1.468E + 5 

7*>. 


3.118E-3 


.V 


4.895E + 4 


The shear stresses on the 30 slip planes are calculated using Eqs. (15-16) as 
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5.995E +4 
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r 16 




T 25 
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- 1.615E + 5 





t 5 


- I.199E + 5 


r 17 




r 26 


-6.922E + 4 







1.154E + 5 





T 6 


-1.599E + 5 


t 18 




r 27 


8.652E + 4 

< 

► = < 


s 4 


> ~ < 

4.615E + 4 

% < 


► = < 

• 

T 1 


-5.995E+4 


r 19 




t 28 


-1.73E + 4 
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-3.996E+4 


T 20 




t 29 


1.557E + 5 







- 9.229E + 4 
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-9.99 IE +4 


f 21 
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- 5.191E + 4 







1.154E+4 


k J 
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( 22 ) 


The engineering shear strains on the 30 slip planes are 


W 1 


- 9. USE - 3 


\r li ] 


-0.014 
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0.017 


y 14 


-1.364E-3 
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7.362E - 3 


y 15 


0.015 



r 4 


8.452E - 3 


y 16 


-0.018 


y*' 


8.8 18E - 3 ' 

f 


-0.011 


r 17 


0.016 


yl & 


- 4.409E - 3 

r 6 


-0.02 


y™ 


1.575E-3 


r 27 


5.51 IE -3 




< 


> — 1 


< 



• 

f 


- 2.742E - 4 


f 9 


0.014 


y 2S 


- 1.102E - 3 

f 


-0.012 


yl 0 


- 7.243E - 3 


y 29 


9.92E - 3 

f 


-0.012 


y 21 


- 6.768E - 3 


j 30 . 


- 3.307E - 3 

y™ 


9.45 IE - 3 


y 22 


-1.364E-3 



y u 


5.907E - 3 


y* 


- 7.502E - 3 





-3.544E-3, 


J 24 . 


8.867E - 3 , 




( 23 ) 



Nagaraj K. Arakere and Gregory R. Swanson 


Page 9 


The normal stresses and strains on the principal and secondary octahedral planes are 

293,700] [ef] [ 0.007185 ' 

130,500 0.001989 

L ^ 2 = J ^ (24) 

32,630 e" -0.001128 

0 el 1-0.002167 

The normal stresses and strains on the cube slip planes are simply the normal stresses and strains 
in the material coordinate system along (100), (010), and (001) axes. This procedure is used to 
compute normal and shear stresses and strains in the material coordinate system, for uniaxial test 
specimens loaded in strain control, in different orientations. 

Strain controlled LCF tests conducted at 1200 F in air for PWA1480/1493 uniaxial smooth 
specimens, for four different orientations, is shown in Table 3. The four specimen orientations 
are <001> (5 data points), <11 1> (7 data points), <213> (4 data points), and <01 1> (3 data 
points), for a total of 19 data points. Figure 5 shows the plot of strain range vs. Cycles to failure. 
A wide scatter is observed in the data with poor correlation for a power law Fit. The first step 
towards applying the failure criteria discussed earlier is to compute the shear and normal stresses 
and strains on all the 30 slip systems, for each data point, for maximum and minimum test strain 
values, as outlined in the example problem. The maximum shear stress and strain for each data 
point, for min and max test strain values, is selected from the 30 values corresponding to the 30 
slip systems. The maximum normal stress and strain value on the planes where the shear stress is 
maximum is also noted. These values are tabulated in Table 4. Both the maximum shear stress 
and maximum shear strain occur on the same slip system, for the 4 different configurations 
examined. For the <00 1> and <01 1> configurations the max shear stress and strain occur on the 
secondary slip system (x 14 , y 14 and x 15 , y 15 respectively). For the <11 1> and <213> configurations 
max shear stress and strain occur on the cube slip system (x 25 , y 25 and x 29 , y 29 respectively). Using 
Table 4 the composite failure parameters highlighted in Eqs. (1-4) can be calculated and plotted 
as a function of cycles to failure. In addition to the four failure parameters discussed, some other 
composite parameters are also plotted as a function of cycles to failure (N). 

Figures (6-9) show that the four parameters based on polycrystalline fatigue failure parameters 
do not correlate well with the test data. The application of these parameters for single crystal 
material is somewhat different since they are evaluated on the slip systems that are thought to be 
the critical planes. The parameters that collapse the failure data well and give the best correlation 

with a power law fit are the maximum shear stress amplitude [Az^ ax ] shown in Fig. 10, the 
composite parameter [(Aw)(A)(nax /2)] shown in Fig. 11, and the max principal shear stress 

amplitude (Tresca theory) shown in Fig. 13. The parameter Az^, a x is appealing to use for its 
simplicity, and its power law curve fit is shown in Eq. (25). It must be remembered that these 
curve fits are only valid for 1200F. 

AT m » = 397,758 N ° 1598 



(25) 
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The correlation for [At^x] would be better if some of the high stress data points are corrected 
for inelastic effects. Since the deformation mechanisms in single crystals are controlled by the 
propagation of dislocations driven by shear the [At^x ] might indeed be a good fatigue failure 
parameter to use. This parameter must be verified for a wider range of R-values and specimen 
orientations, and also at different temperatures and environmental conditions. Equation 25 will 
be used to calculate fatigue life at a critical blade tip location for the SSME turbine blade. 


5.0 Application of Fatigue Failure Criteria to Finite Element Stress Analysis Results of 
Single Crystal Nickel Turbine Blades 

Turbine blades used in the advanced high-pressure fuel turbopump (AHPFTP) are fabricated 
from single crystal nickel PWA1480/1493 material. Many of these blades have failed during 
operation due to the initiation and propagation of fatigue cracks from an area of high 
concentrated stress at the blade tip leading edge. Inspection of blades from other units in the test 
program revealed the presence of similar cracks in the turbine blades. During the course of the 
investigation an interesting development was brought to light. When the size of the fatigue 
cracks for the population of blades inspected was compared with the secondary crystallographic 
orientation (3 a definite relationship was apparent as shown in Fig. 14 2> 14 . Secondary orientation 
does appear to have some influence over whether a crack will initiate and arrest or continue to 
grow until failure of the blade airfoil occurs. Figure 14 reveals that for (3 = 45+/- 15 deg tip 
cracks arrested after some growth or did not initiate at all. This suggests that perhaps there are 
preferential [3 orientations for which crack growth is minimized at the blade tip. 

In an attempt to understand the effect of crystal orientation on blade stress response a three- 
dimensional FE model capable of accounting for primary and secondary crystal orientation 
variation was constructed. The Alternate High Pressure Fuel Turbo Pump (HPFTP/AT) first 
stage blade ANSYS finite element model (FEM) was cut from a large 3D cyclic symmetry 
model that includes the first and second stage blades and retainers, interstage spacer, disk and 
shaft, and the disk covers (Fig. 15). The blade dampers are represented with forces applied to the 
blade platforms at the damper contact locations. The models are geometrically nonlinear due to 
the contact surfaces between the separate components. The element type used for the blade 
material is the ANSYS SOLID45, an 8-noded 3D solid isoparametric element. Anisotropic 
material properties are allowed with this element type. ANSYS aligns the material coordinate 
system with the element coordinate system. 

The effect of crystal orientation on blade stress response was studied by running 297 separate FE 
models to cover the complete range of primary and secondary crystal orientation variation. To 
generate the 297 material coordinate systems used for this study local coordinate systems were 
generated and the element coordinate systems aligned with them 15 . The material coordinate 
system is referenced to the blade casting coordinate system, shown if Fig. 16. The relative 
orientation of the primary axis of the material coordinate system to the casting coordinate system 
is shown in Fig. 17. Two angles, A and y, locate the primary material axis relative to the casting 
axis, the third angle, (3, is the clocking of the secondary material axis about the primary material 
axis. Figure 18 and Table 5 show the distribution of the 297 different material coordinate 
systems within the allowed 15 degree maximum deviation from the casting axis. The secondary 
repeats after 90 degrees, so only 0 to 80 degrees needs to be modeled. 
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The load conditions represent full power mainstage operation of the Space Shuttle Main Engine, 
referred to as 109% RPL SL (Rated Power Level Service Life). The shaft speed is 37,355 rpm, 
the airfoil temperature is approximately 1200 F, forces representing the blade damper radial sling 
load are applied to the blade platform, and aerodynamic pressure loads are applied to the blade 
surfaces and internal core. 

Postprocessing of the 297 finite element results files presented a fairly difficult challenge, and 
represents a considerable amount of effort. Two FORTRAN programs were employed for the 
postprocessing work. The first selects the element results from the coded binary output files and 
places them into ASCII text files. The second program processes the ASCII files to calculate 
averaged nodal results, the resolved shear stresses and strains and the normal stresses and strains 
on the 30 slip systems, in the single crystal material coordinate system. It then calculates the 
parameters chosen for study and sorts them based on user set criteria. 

The connection between the blade and disk are modeled with ANSYS COMBIN40 elements. 
These elements have one degree of freedom at each node. The nodal motion in that degree of 
freedom sets the separation or contact of these elements only. This element does not have the 
capability for friction tangent to the contact surface. For this model the nodal coordinate systems 
on the contacting surface of the blade firtree attachment were rotated so that one axis is normal 
to the surface. This is the degree of freedom used in the COMBIN40 element. The nodal 
coordinate systems on the disk contact surfaces were similarly oriented. An interesting feature 
of the ADAPCO model is that the blade is next to a cyclic symmetry section of the disk (a 1 of 
50 piece) so that only the pressure side of the blade attachment contact surface nodes are nearly 
coincident to the disk. The suction side of the blade is clocked 7.2 degrees about the shaft from 
the mating surface on the disk. The blade and disk nodal coordinate systems for the suction side 
attachment are 7.2 degrees out of parallel to each other to account for this. Since the 
COMBIN40 element only acts on the one degree of freedom normal to the contact surfaces the 
7.2 degree offset in physical location and orientation is properly accounted for. To run the blade 
model separate from the global model the nodal displacements of the disk nodes attached to the 
COMBIN40 elements were taken from a run of the global model and used as enforced 
displacements for what would become free ends of the contact elements. 

Figure 19 shows representative Von Mises stress distribution plot for the turbine blade in the 
attachment region. The crack location and orientation at the critical blade tip location is shown in 
Fig. 14. 


6 Effect of Secondary Crystal Orientation on Blade Tip Stress Response 

Variation of secondary crystal orientation on stress response at the blade tip critical point prone 
to cracking (tip point on inside radius) was examined by analyzing the results from the 297 FE 
model runs. The FE node at the critical point was isolated and critical failure parameter value 
(Ax inax ) computed on the 30 slip systems. A shaded contour plot of Ax niax was generated as a 
function of primary and secondary orientation, shown in Fig. 20. The contour plot clearly shows 
a minimum value for Ax [n a X for secondary orientation of P = 50 deg and primary orientation 
designated by cases 5 and 20. From Table 5.1 we see that case 5 corresponds to a primary 
orientation of A = 0 deg and y = 7.5 deg. Case 20 corresponds to a primary orientation of A = 
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5.74 deg and y = 13.86 deg. Using the fatigue life equation based on the AXmax curve fit of LCF 
test data, Eq. (25), we can obtain a contour plot of dimensionless life at the critical point as a 
function of primary and secondary orientation, as shown in Figure 21. The maximum life is 
again obtained for P = 50 deg, and A = 0 deg and y = 7.5 deg, and A = 5.74 deg and y = 13.86 
deg. The optimum value of secondary orientation P = 50 deg. corresponds very closely to the 
optimum value of P indicated in Fig. 14. This demonstrates that control of secondary and 
primary crystallographic orientation has the potential to significantly increase a component’s 
resistance to fatigue crack growth without adding additional weight or cost. 


7.0 CONCLUSIONS 

Fatigue failure in PWA1480/1493, a single crystal nickel base turbine blade superalloy, is 
investigated using a combination of experimental LCF fatigue data and 3D FE modeling of 
HPFTP/AT SSME turbine blades. Several failure criteria, based on the normal and shear stresses 
and strains on the 24 octahedral and 6 cube slip systems for a FCC crystal, are evaluated for 
strain controlled uniaxial LCF data (1200 F in air). The maximum shear stress amplitude [Armax] 
on the 30 slip systems was found to be an effective fatigue failure criterion, based on the curve 
fit between Ar™* and cycles to failure. Since deformation mechanisms in single crystals are 
controlled by the propagation of dislocations driven by shear, Ar^x might indeed be a good 
fatigue failure parameter to use. However, this parameter must be verified for a wider range of 
R-values and specimen orientations, and also at different temperatures and environmental 
conditions. 

Investigation of leading edge tip cracks in operational SSME turbine blades had revealed that 
secondary crystal orientation appeared to influence whether a crack initiated and arrested or 
continued to grow until failure of the blade airfoil. The turbine blade was modeled using 3D FEA 
that is capable of accounting for material orthotrophy and variation in primary and secondary 
crystal orientation. Effects of variation in crystal orientation on blade stress response were 
studied based on 297 FE model runs. Fatigue life at the critical locations in blade was computed 
using FE stress results and failure criterion developed. Detailed analysis of the results revealed 
that secondary crystal orientation had a pronounced effect on fatigue life. The optimum value of 
secondary orientation (5 = 50° computed corresponds very closely to the optimum value of (3 
indicated in the failed population of blades. Control of secondary and primary crystallographic 
orientation has the potential to significantly increase a component’s resistance to fatigue crack 
growth without adding additional weight or cost. 
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Table 3 Strain controlled LCF test data for PWA1493 at 1200 F for 4 specimen orientations 





































= Max shear strain of 30 slip systems for max specimen test strain value 
= Max shear strain of 30 slip systems for min specimen test strain value 
= Max shear stress of 30 slip systems for max specimen test strain value 
= Max shear stress of 30 slip systems for min specimen test strain value 
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Case 

Delta 

Gamma 

Beta 

0 

0.00 

0.00 

0, 1 0,20,30,40,50,60,70,80 

1 

7.50 

0.00 

0, 1 0,20,30,40,50,60,70,80 

2 

6.93 

2.87 

0,10,20,30,40,50,60,70,80 

3 

5.30 

5.30 

0,10,20,30,40,50,60,70,80 

4 

2.87 

6.93 

0,10,20,30,40,50,60,70,80 

5 

0.00 

7.50 

0,10,20,30,40,50,60,70,80 

6 

-2.87 

6.93 

0,10,20,30,40,50,60,70,80 

7 

-5.30 

5.30 

0,10,20,30,40,50,60,70,80 

8 

-6.93 

2.87 

0,10,20,30,40,50,60,70,80 

9 

-7.50 

0.00 

0,10,20,30,40,50,60,70,80 

10 

-6.93 

-2.87 

0,10,20,30,40,50,60,70,80 

11 

-5.30 

-5.30 

0,10,20,30,40,50,60,70,80 

12 

-2.87 

-6.93 

0,10,20,30,40,50,60,70,80 

13 

0.00 

-7.50 

0,10,20,30,40,50.60,70,80 

14 

2.87 

-6.93 

0,10,20,30,40,50,60,70,80 

15 

5.30 

-5.30 

0,10,20,30,40,50,60,70,80 

16 

6.93 

-2.87 

0,10,20,30,40,50,60,70,80 

17 

15.00 

0.00 

0,10,20,30,40,50,60,70,80 

18 

13.86 

5.74 

0,10,20,30,40,50,60,70,80 

19 

10.61 

10.61 

0,10,20,30,40,50,60,70,80 

20 

5.74 

13.86 

0, 1 0,20,30,40,50,60,70,80 

21 

0.00 

15.00 

0,10,20,30,40,50,60,70,80 

22 

-5.74 

13.86 

0,10,20,30,40,50,60,70,80 

23 

-10.61 

10.61 

0, 1 0,20,30,40,50,60,70,80 

24 

-13.86 

5.74 

0,10,20,30,40,50,60,70,80 

25 

-15.00 

0.00 

0, 1 0,20,30,40,50,60,70,80 

26 

-13.86 

-5.74 

0,10,20,30,40,50,60,70,80 

27 

-10.61 

-10.61 

0, 1 0,20,30,40,50,60,70,80 

28 

-5.74 

-13.86 

0,10,20,30,40,50,60,70,80 

29 

0.00 

-15.00 

0,10,20,30,40,50,60,70,80 

1 30 

5.74 

-13.86 

0,1 0,20,30,40,50,60,70,80 

31 

10.61 

-10.61 
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32 

i 13.86 

-5.74 

0,10,20,30,40,50,60,70,80 


Table 5 33 primary axis cases with 9 secondary cases each, for a total of 
297 material orientations. 





































































































Fig. 1 Convention for Defining Crystal Orientation in Turbine Blades 2 















Strain Amplitude (Uniaxial LCF Data) 
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Power Law Curve Fit (R A 2 = 0.469): As = 0,0238 N 
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Fig. 5 Strain range Vs. Cycles to Failure for LCF test data (PWA1493 at 1200F) 
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Power Law Curve Fit (R A 2 = 0.130 ): [ft*,* + £ n ] = 0.0249 N' 0 ' 773 



Fig. 6 [%iax "F &n ] (Eq. 1) Vs. N 


Delta (Gamma+Eps)/2 + SigmaMean/E = Eq. (3.2) 


Power Law Curve Fit (R A 2 = 0.391 ): 
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Power Law Curve Fit (R A 2 = 0.674 ): Ax = 397,758 N 0 1598 



Fig. 10 Shear Stress Amplitude [Av ax ] Vs. N 
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Power Law Curve Fit (R A 2 = 0.549 ): T max * Ay/2 = 4,661 N 
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Figure 14 Secondary Crystallographic Orientation, |3, Vs Crack Depth for the SSME 
AigS^TP 1st Stage Turbine Blade 2 
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Figure 18 33 primary axis cases with 9 secondary cases each, for a total of 297 material orientations. 
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Figure 16. First Stage Blade FE Model and Casting 
Coordinate System Coordinate System 




Figure 19 Representative von Mises stress distribution results in the blade attachment region 
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Fig. 20 Maximum Shear Stress Amplitude (Alma*, ksi) Contour at the Blade Tip Critical Point 
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